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Abstract 

Metallothioneins (MTs) are excellent candidate genes for Inflammatory Bowel Disease (IBD) and have previously 
been shown to have altered expression in both animal and human studies of IBD. This is the first study to examine 
genetic variants within the MT genes and aims to determine whether such genetic variants have an important role 
in this disease. 28 tag SNPs in genes MT1 (subtypes A, B, E, F, G, H, M, X), MT2, MT3 and MT4 were selected for 
genotyping in a well-characterized New Zealand dataset consisting of 406 patients with Crohn's Disease and 638 
controls. We did not find any evidence of association for MT genetic variation with CD. The lack of association 
indicates that genetic variants in the MT genes do not play a significant role in predisposing to CD in the New 
Zealand population. 

Keywords: Crohn's disease, Metallothioneins, genetic association 



Background 

Inflammatory bowel disease (IBD) is a disorder charac- 
terised by chronic, relapsing inflammation of the gastro- 
intestinal tract and is characterized by the presence of 
an increased level of reactive oxygen species in the 
mucosal intestinal tissue as well as oxidative DNA and 
protein damage, defective host-microbe interactions, 
immune cell infiltration, and disturbed T cell apoptosis. 
Metallothioneins (MTs) are able to affect all of these 
processes, making them good candidates for IBD. MTs 
are ubiquitous metal-binding proteins that have been 
highly conserved throughout evolution and are rapidly 
upregulated in response to an inflammatory stimulus. 
There are four main isoforms expressed in humans: 
MT1 (subtypes A, B, E, F, G, H, M, X), MT2, MT3 and 
MT4. These genes cluster together on a single locus on 
chromosome 16 (16ql3) [1]. 

Support for the hypothesis of MTs playing an 
important role in the pathogenesis of IBD has come 
from various reports demonstrating altered MT 
expression in IBD. However these reports have been 
contradictory (for a summary see Waeytens et al., 
2009 [2]) with some studies reporting MT 
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upregulation in IBD [3,4] and others demonstrating a 
downregulation [5-11]. The inconsistencies between 
the various studies may be explained by differences in 
patients samples such as age, medication, disease 
activity, zinc status, as well as from where the tissue 
was sampled (ileum or colon). 

Whilst there have been several studies investigating 
expression of MTs in IBD, there have been no studies 
published to date investigating genetic variants in MTs 
and IBD. The aim of this study is to identify SNPs 
across these genes and to genotype these in a well-char- 
acterised Caucasian New Zealand IBD dataset and to 
examine the results for evidence of genetic association 
between any of the MT genes and IBD. 

Methods 

Samples 

A total of 1044 subjects from New Zealand were 
included in the study: 406 CD patients and 638 age and 
sex matched controls. All participants self-reported Eur- 
opean ancestry. 

Clinical records were analysed to confirm diagnosis, 
and IBD status was defined using standard diagnostic 
criteria [12]. Cases were phenotyped according to the 
Montreal Classification systems. Clinical characteristics 
of the CD patients are shown in table 1. 
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Table 1 Summary of clinical data of CD patients. 
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Participants consented to collection of peripheral 
blood or a buccal swab for DNA extraction and geno- 
typing, and DNA was extracted from the blood/buccal 
samples using Qiagens DNA extraction kit and follow- 
ing the manufacturer's instructions. 

The study was conducted under ethical protocol 
MEC/04/12/011, authorised through the New Zealand 
Multi- Region Human Ethics Committee. All study sub- 
jects gave informed consent. 

Snp selection 

Tag SNPs in genes MT1 (subtypes A, B, E, F, G, H, M, 
X), MT2, MT3 and MT4 were selected using Hapmap 
release 28, and the tagger functionality within Haploview 
with pairwise tagging to identify SNPs using an r 2 > 0.8 
and a minor allele frequency > 5%. As a result 28 tag 
SNPs were selected for genotyping. 

Genotyping 

Genotyping was performed with the MassARRAY and 
iPlex systems of the Sequenom genotyping platform 
(Sequenom, San Diego, CA), which uses the MALDI- 
TOF primer extension assay [13,14], according to manu- 
facturers' recommendations. 

Assays were optimized in 24 samples consisting of 20 
reference Centre d'Etude du Polymorphisme Humain 
(CEPH) samples and 4 blanks. 

All sample plates contained cases, controls, blanks, 
CEPH and duplicate samples. Quality control measures 
included independent double genotyping and, where 



available, comparison of our CEPH genotypes to those 
in the Hapmap database http://www.hapmap.org. 

Statistical analysis 

SNPs were tested for deviation from HWE in both cases 
and controls using a chi-square goodness-of-fit test. To 
determine if there were differences between cases and 
controls, genotype and allele frequencies for each SNP 
were analyzed by 2 x 3 and 2x2 Chi-square tables 
respectively 

Genotype and phenotype associations were assessed by 
comparing allele frequencies between controls and 
patient subgroups defined using the clinical characteris- 
tics. These analyses were carried out using R (R: A lan- 
guage and environment for statistical computing, R 
Foundation for Statistical Computing, Vienna, Austria. 
ISBN 3-900051-07-0, URL http://www.R-project.org/) 
and SAS (V9.1 SAS Institute., Cary, NC, USA). 

To determine linkage disequilibrium (LD) between 
SNPs and to define haplotype blocks, we uploaded our 
data into Haploview [15]. Haplotype blocks were defined 
using the default algorithm which uses confidence inter- 
vals [16]. Haplotype analysis was carried using HAPLO. 
SCORE in R to test for association of these haplotypes 
with CD. 

For all analyses we considered a p values less than 
0.05 to indicate statistical significance. 

The false discovery rate (FDR) was used to correct for 
multiple testing [17,18]. 

Results and discussion 

Three SNPs did not fit into the Sequenom multiplexes: 
rsl827210 (MT1M), rsll076161 (MT1A), rs2298846 
(MT1G). Two SNPs failed to be successfully geno- 
typed: rsl827208 (MT1M) and rsl2448654 (MT1G). 
One SNP was out of HWE in the control samples and 
so was excluded from further analysis: rsl2051311 
(MT1B). Thus there were 22 SNPs with genotype data 
for analysis. Genotype and allele counts/frequencies 
and p-values are shown in table 2. One SNP rs4784708 
(MT1H) was found to be associated with CD (genoty- 
pic p = 0.0197, allelic p = 0.016) but after applying 
multiple testing correction using FDR the result is no 
longer statistically significant. Thus we can conclude 
that individual SNPs in the MT genes are not asso- 
ciated with CD. 

Calculations of statistical power using PS 2.1.31 [19] 
show that using our case-control sample set although we 
have good power to detect an odds ratio of 1.3 for com- 
mon SNPs (for example for a SNP with MAF 0.42 (e.g., 
rs 11644094) we have 86.5% power) the power is consid- 
erably lower for SNPs with low MAF. For example for a 
SNP with MAF of 0.02 (e.g., rs9936741) we have 16.1% 
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Table 2 Genotype and allele counts (and frequencies) in CD patients and in controls 
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Table 2 Genotype and allele counts (and frequencies) in CD patients and in controls (Continued) 
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power. Thus, an association of metallothionein genes 
with CD cannot be entirely ruled out as it may simply be 
that we do not have the power in our study to detect 
such an association. A replication study in a larger sam- 
ple set is encouraged. 

There are 5 haplotype blocks in this region (Figure 1) 
but examining the different haplotypes we did not find 
any associations with CD (data not shown but available 
on request). 

We also undertook a phenotype analysis and after 
multiple testing corrections none of the SNPs were 
associated with any of the different clinical characteris- 
tics of CD (data not shown but available on request). 

This is the first study to investigate genetic association 
with MT SNPs and IBD. In summary although MT 
genes are good candidates for CD, the present study 
failed to detect any significant associations. The lack of 
association indicates that genetic variants in the MT 
genes do not play a significant role in predisposing to 
CD in the New Zealand population. 
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